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Abstract-Genetic algorithm (GA) is an intelligent method based 
on natural selection for global numerical optimization. To make 
use of the great advantages of GA, this paper focuses on applying 
GA to multi-dimensional bin-packing problem (MBP). In multi- 
dimensional bin-packing problem, the primary goal is putting 
grouping items into appropriate bins with the same dimensions to 
make full use of the bins and minimize the number of bins used. 
In other words, the optimal solution requires that most bins are 
filled in maximum degree. By designing proper fitness functions, 
GA is successfully applied to MBP. The experimental results 
reveal that the proposed GA method solves the MBP efficiently. 
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I. INTRODUCTION 

During the past decade, bin-packing (BP) problem [1, 2, 3] 
has become a research hotspot in combinatorial optimization. A 
lot of researchers have devoted to the theoretical developments, 
thus many heuristic algorithms are explored for the bin-packing 
problem. However, the topics are far from ending, there are no 
mature theory and exact method up to now. Under this 
circumstance, researchers mainly use approximate heuristic 
algorithm to deal with BP problem. Johnson analyzed the 
performance of the First-Fit (FF), Best-Fit (BF), First-Fit 
Decreasing (FFD) and Best-Fit Decreasing (BFD) heuristics in 
[2], Bhatia proposed a Better-Fit Heuristic Method for One- 
dimension BP Problem in [4]. In recent years, more work has 
been done to develop the algorithms for two-dimension and 
three-dimension BP problem [5, 6], Besides that, intelligent 
methods are applied to BP problem to optimal the above 
methods, such as, ACO, PSO. 

Taking the process of rotation in real packing situation into 
consideration, Genetic algorithm (GA) is the most common 
used one in MBP which can simulate the process of rotation by 
designing suitable chromosome and fitness function. GA [7] is 
an efficient method for global optimization. So, it is suitable to 
generalize GA method to BP problem. There are many 
applications of GA for BP problem. A Grouping Genetic 
Algorithm was developed by Falkenauer in [6]. And in [8], 
Falkenauer developed a new representation and operators of 
GA for BP problem. 

In this paper, due to the sort factor in BFD is special, an 
improved BFD (IBFD) will be introduced. The aim of this 
paper is to develop an improved heuristic algorithm based on 
GA (MBPGA) which will be introduced in details later. This 
algorithm discusses the way of solving the MBP with 
dimensions bigger than 3 (d>3). The goal of this algorithm is to 
improve IBFD to seek the optimal solution that the sum of 
every dimension of the grouping items in each bin is less than 
or equal to the capacity of bin. In this paper, the predefined 
bins and grouping items are all rotatable. 



The rest of the paper is organized as follows. Some existed 
heuristic algorithms together with definitions for MBP (d>3) 
are presented in Section 2. Then, a comparison is presented 
between IBFD method and MBPGA which lead out the way 
defining fitness functions for the MBP in Section 3. Within 
Section 4, a practical application for MBP (d>3) is given. 
Finally, conclusions and future work are provided in Section 5. 

II. HEURISTIC ALGORITHM AND DEFINITIONS 
FOR MBP 

In this section, some existed heuristic algorithms for BP 
problem will be introduced and some definitions for MBP are 
stated. 

A. Some Popular Heuristic Algorithm 

Various heuristic on-line algorithms as well as off-line 
methods have been applied to BP problem by now [9], 
Generally, on-line algorithms pack the items to a bin in 
sequence of arrival while off-line methods have all items 
available before packing. Three on-line (NF, FF, BF) together 
with two off-line methods are introduced as follows. 

1) On-line algorithms 

• Next-Fit (NF): NF keeps a bin active permanently, it 
puts the item into the bin if it fits in the bin; otherwise the 
heuristic uses a new bin for the item which is marked as the 
active bin. 

• First Fit (FF): FF checks all the non-empty bins and 
packs the item into the first bin it fits; otherwise FF packs the 
item in a new bin. 

• Best Fit (BF): BF packs the item into a bin that is 
filled to maximum degree but still having enough vacant space 
to pack the item; otherwise BF packs the item in a new bin. 

2) Off-line algorithms 

• First-Fit Decreasing (FFD): FFD packs the items 
according to the FF after sorting the list of items in non- 
increasing order of sizes. 

• Best-Fit Decreasing (BFD): BFD packs the items 
according to the BF after sorting the list of items in non- 
increasing order of sizes. 

Among the above methods, BFD is the best one with most 
efficiently packing effect. This is due to its sorting method to 
packing the biggest item into the smallest bin which has 
enough vacant space, which makes it possible to achieve a goal 
of making the most use of available bins. 

B. Concepts and Definitions for MBP 

In this part, some definitions and the Improved BFD (IBFD) 
are introduced. 
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Definition 1 

A sign function, sign(x), is defined by 

' 1 ,x>0 
sign(x) 



{ ,x = 
l,x<0 



(1) 



Definition 2 



Let scalar d be the dimensions of the items and bins, m be 
the number of items and n be the number of bins. 

Let ITEMS shown in Eq. (2) be a matrix which consists of 
it tj (1< i < d,l< j < m) . The column number of ITEMS 

denotes the number of each item and the row number denotes 
the number of dimension corresponding to each item. For 

example, it tj is the element in row i and column j which 
denotes the size of i-th dimension in item j. 



ITEMS 



... it ±] . 


■ 


... it 2] . 


•• a 2m 


... it v . 


•■ K 


••• Kj ■ 


■■ K m , 



(2) 



Note that, IT 1 is used to represent the column vector i i.e. 
TT i =(it il ,it 2l ,.Jt dl ) T . 

Let BINS shown in Eq. (3) be a matrix which consists of 
b tj (l<i<d,l< j <n) . The column number of BINS 

denotes the number of each bin and the row number denotes 
the number of dimension corresponding to each bin. For 

example, £>.. is the element in row i and column j which 
denotes the size of i-th dimension in bin j. 



BINS 






\b dl 



(3) 



dn J 



Note that, B J is used to represent the column vector j i.e. 
B j =(b 1 .,b 2j ,...b dj ) T 

Definition 3 

In this paper, MBP (d>3) only compares the size of every 
dimension in the items and bins. If the sum of every dimension 
of the grouping items in each bin is less than or equal to the 
capacity of bin, IBFD packs the item into the bin, In other 

words, every element in IT must be less than or equal to the 
corresponding element in B 1 . For example, if 
IT =(0.2, 0.5, 0.3) r and B J =(0.5,0.6,0.8) r , it is 
obvious that every element in IT' is less than the 

corresponding element in B 1 , then IBFD packs the item i into 
the bin j. 



Let STATES shown in Eq. (4) be a matrix which consists of 
S y (1< i < m,l< j < n) . The column number of STATES 
denotes the number of each bin and the row number denotes 
the number of each bin. For example, S. is the element in row 
i and column j. If IBFD packs the item i into the bin j, then 



STATES 



a ll a 12 ' ' ■ J lj 

s 21 s 22 ... S 2 j 

S n S [2 ... Sy 


• s l„ 

• S 2„ 



(4) 



(5) 



Definition 4 

In IBFD, the sort factor is defined by 

sort factor = the sum of every dimmensionx 
theproductof every dimmension 

For example, if IT 1 = (0.2, 0.5, 0.3) r , then 

sort factor = (0.2 + 0.5 + 0.3) x(0.2 x 0.5 x 0.3) 
= 0.03 

IBFD packs the items according to the BF after sorting the 
list of items in non-increasing order and the list of bins in non- 
decreasing order of each sort factor. 

III. DESIGN OF MBPGA FOR MBP 

In this section, MBPGA is applied to MBP and makes a 
comparison with the traditional IBFD. 

A. Packing Items Based on IBFD 

Let's assume that ITEMS shown in Eq. (2) and BINS shown 
in Eq. (3) are already sorted by sort rules which are introduced 

above, then IT is the largest items and B is the smallest bin 

in each list. First, IBFD packs IT into first bin that having 

enough vacant space to pack it, such as B'. After packing, 

B 1 = (b n -it u ,b 2t -it 21 , ...,b di -it dl ) T , then IBFD sorts the 

bins again and packs the next item. Sometimes, the next item 

(IT ) cannot packs into B' for at least one dimension in IT 

is bigger than the corresponding dimension in B . However, 

considering the item is rotatable, we can rotate IT , then IT 

is packed into B' . Unfortunately, in IBFD, there doesn't exist 
a rotate process, so some bins are not filled to the maximum 
degree. To solve this problem, MBPGA is developed in the 
next part. 

B. Solving MBP by MBPGA 

MBPGA is a heuristic algorithm that improves IBFD based 
onGA. 

Definition 5 

Let Q m denote all the orders of items' with m dimensions 
for packing. For example, Q 4 = (1, 2, 3, 4) is an initial 
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sequence and Q 4 =(2,3,1,4) denotes the sequence 
corresponding to the bin after the process of rotation. 

Let IT i denotes the i-th member in Q n 
Q 4 =(2,3,1,4) .then IT 2 =3. 



1) Procedures ofMBPGA 

To apply the MBPGA to this case, the procedures are as 
follows [10]. 

a) Initialize a Population. 

Each order of item's dimensions is a member of the 
population. In MBPGA, the order of item's dimensions is 
called a chromosome which is coded in real number in this 
paper. To make it clear for understanding, taking d-A to 
demonstrate a chromosome. {3,1,4,2} is a chromosome which 
represents the sequence corresponding to the bin. 

In this case, because the initial population are generated 
randomly, there may be some invalid chromosomes such as 
one dimension is selected for more than one time, for example, 
{2,1,2,1}, {1,3,3,3}, and {4,4,4,4}. 

b) Define the Fitness Function. 

How to define fitness function to exclude these invalid 
chromosomes and, at the same time, to express the fitness 
value of every member exactly, is a crucial process. This 
process will be discussed in details later. 

c) Calculate the Fitness Value. 

d) Selection. 

Selecting a new population based on the fitness value above. 

e) Crossover. 

Single point crossover is applied to this process, such as, 
two chromosomes {1,2,5,3} and {4,2,3,1}, after crossover, 
these two chromosomes may change to {1,2,3,1} and {4,2,5,3}. 

f) Mutation. 

In this process, a random probability is occurred to 
determine whether there exist changes in the chromosome. 

g) Repeat the Processes c to f until the Best Packing 
Solution Is Worked out. 

2) Fitness Function for MBPGA 

The fitness function of MBPGA takes the role of the natural 
environment to screening the inferior individuals and retaining 
the fitted one, so it is a key part in MBPGA. 

In this case, fitness function is defined as 

fitness _value = ]~J sign(\B' -IT \+(B' -IT'))- 



is[l,m] 
Ml,n] 



(6) 



signill \IT l -IT j \)-(V £ B'-IT') 



i,j<m m 



ie[l,m] 

Ml."] 



In Eq. (6), Y[ sign(\ B 1 - IT' \+(B ] - IT')) 



is 



i«s[l,m] 
je[l,n] 



designed to screening the situation that, at least in one 
dimension, the item is bigger than the bin. In that situation, 

there exists B J < IT' , then sign(\ B> -IT 1 \ +(B j -IT')) = 

as sign(x) is expressed in Eq. (1). The second component is 



constructed to get rid of the invalid chromosomes. If there 
exists one or more same dimension numbers in Q m , then 



such as signCjJ I T, ~ ^ |) = ■ Finall y> £ b j _ IT < 



is the 



i,J6S!„ 

size of the bin after packing, so 



is[l,m] 



a/Z 



IT") 



denotes 



etl.m] 

Ml,»] 



the degree of utilization ratio of the bin. If one of the first two 
components equal zero, then the fitness value will be zero, in 
other words, this chromosome will not be remained in the next 
generation for its fitness value is small. 

IV. A PRACTICAL APPLICATION 

In this section, a practical application of MBP (d=4, m=2, 
n=2) is taken as an example to demonstrate the way MBPGA 
works. Moreover, a comparison is made between the 
performances of these two methods by applying them to MBP 
(d=4, m=6, n=4). 

ITEMS and BINS are initialized as (7) and (8). 

f0.15 0.2^ 
0.35 0.3 



ITEMS 



BINS 



0.6 


0.7 


0.9 


0.8 


A 1.8 


1.2 


1.1 


1.1 


1.0 


1.05 


,1.2 


1.8 



(7) 



(8) 



A. Using IBFD 

In this part, IBFD is applied to MBP as follows: 

1) Calculate the sort factors of the ITEMS and BINS 
factor of IT X = (0.15 + 0.35 + 0.6 + 0.9)*(0.15*0.35*0.6* 0.9) = 0.0567 
factor of IT 2 = (0.2 + 0.3 + 0.7 + 0.8)*(0.2*0.3*0.7*0.8) = 0.0672 

factor of B l = (1+1.1+1.2 +1,8)*(1*1.1*1,2*1,8) =12.1176 
factor of B 2 = (1.05 + 1. 1 + 1. 2 +1.8)*(1.05*1. 1*1. 2*1.8) =12.848 

2) Sorting the List of Items in Non-Increasing Order and 
the List of Bins in Non-Decreasing Order 

After sorting, the list of items and bins are shown in (9) and 
(10). 



ITEMS 



BINS = 



(0.2 


0.15^ 


0.3 


0.35 


0.7 


0.6 


v 0.8 


0.9 , 


'1.8 


1.2^ 


1.1 


1.1 


1.0 


1.05 


1.2 


1.8 J 



(9) 



(10) 



3) Packing 

Firstly, IBFD packs IT into B , then if B is still having 
enough vacant space to pack IT , IT will be packed into 
B ; otherwise IBFD packs IT into B . 
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Now, we pack IT into B . After packing, BINS changes 
to (11). 



BINS 



(1.6 


1,2 1 


0.8 


1.1 


0.3 


1.05 


I 0.4 


1.8 J 



(11) 



ITEMS 



It is obvious that B doesn't have enough vacant space to 

pack IT for the third and the forth dimensions of B are 

smaller than IT , in this situation, IT is packed into B . 
However, this solution wastes a lot capacity of both two bins. 

After packing, IT is packed into B and IT is packed 
into B . Two bins are all used in this packing. 

B. Using MBPGA 

To minimize the number of used bins, MBPGA is 
introduced to this case. In MBPGA, the first two steps are the 
same as the one in IBFD. So, here we start from step3. 

Before packing, MBPGA rotates IT . After this process, 
ITEMS is shown in (12). 

^0.8 
0.3 0.35 (12) 

0.2 
v 0.7 

Then IT is packed into B . After packing, BINS changes 
to (13). 

f 1 

0.8 

0.8 
^0.5 

MBPGA will rotate IT before it is packed. After rotation, 
ITEMS is shown in (14). 

^0.8 

0.3 0.35 (14) 

0.2 
x 0.7 

Now, IT is packed into B for B is still having enough 
vacant space to pack IT . 

A conclusion can be drawn from this application that 
MBPGA is an efficient heuristic algorithm for MBP and it is 
more effectively than IBFD in this case. 

C. The Performances of IBFD and MBPGA 

In this part, a practical bin-packing application is 
introduced to illustrate the performances of IBFD and MBPGA. 
Table 1 shows the data of 4 bins, including the size of every 
dimension and the size order of the initial bins. 

TABLE 1 DATA OF BINS 



BINS 



ITEMS 



0.15 A 
0.35 

0.6 
0.9 



1.2^1 

1.1 
1.05 
1.8 



0.9 ^ 

0.35 
0.6 
0.15, 



(13) 



Title 



Size 



Size order (The Min. 
isl) 



Binl 1.0012X1.1888X1.3164X1.6996 

Bin 2 1.2874X1.4390X1.5431X1.6253 

Bin 3 1.5017X1.5761X1.7615X1.7624 

Bin4 1.1232X1.5044X1.6455X1.7477 



The data of 6 items are depicted in Table 2, which also 
consist of the size of every dimension and the size order of the 
given items. 



TABLE 2 DATA OF ITEMS 



Title 



Size 



Size order (The 
Min. is 1) 



Item 1 
Item 2 
Item 3 
Item 4 
Item 5 
Item 6 



0.092148x0.14785x0.19817x0.34726 

0.25678x0.43151x0.67227x0.6944 

0.0097586x0.27939x0.53228x0.94623 

0.024855x0.39268x0.67144x0.90644 

0.056933x0.45032x0.83717x0.9715 

0.58247x0.65004x0.68664x0.71943 



From Table 1 and Table 2, in order to make best use of the 
given bins, it is obvious that the candidate method prefers to 
pack the biggest item (Item 6) into the smallest bin (Bin 1) first, 
if Bin 1 has enough vacant space to pack Item 6. Then, the 
method follows its packing sequence until the smallest item 
(Item 3) is packed. 

Table 3 shows the experiment results of IBFD and MBPGA, 
the performances of these two candidate methods are apparent 
from the results. 

TABLE 3 RESULTS OF TWO METHODS 



Algorithm Bin No. 



Item No. in each 
bin 



Usage rate of each bin (%) 



IBFD 



3,4 

5 

2 
1,6 



99.59 
74.72 
52.82 

91.13 



MBPGA 



1,6 

3,4 

unused 

2,5 



87.85 
89.28 
unused 
95.82 
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IBFD packs all items into the given bins, however, it uses 
up all available bins. Some of the bins are used inefficiently. 
Fortunately, MBPGA is applied to this application by 
designing suitable fitness functions. From Table 3, it is obvious 
that MBPGA uses only 3 bins to pack all the target items. 
Moreover, the biggest bin (Bin 3) is never used which can 
packs other newly arrived items. 

V. CONCLUSIONS AND FUTURE WORK 

In this paper, a novel concept of sort factor and an 
improved BFD (IBFD) are proposed. Then, this paper 
concentrates on designing a new heuristic algorithm to solve 
the MBP (d>3) by improving the process of IBFD. Besides that, 
MBPGA is developed to improve the process of IBFD. 
Moreover, a comparison is made between IBFD and MBPGA 
to analyze the performances of the two algorithms. Finally, 
MBPGA is applied to MBP successfully and efficiently. 
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